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Abstract. This paper deals with the investigation of the solution of an unified 
fractional reaction-diffusion equation associated with the Caputo derivative as 
the time-derivative and Riesz-Feller fractional derivative as the space-derivative. 
The solution is derived by the application of the Laplace and Fourier transforms 
in closed form in terms of the H-function. The results derived are of general 
nature and include the results investigated earlier by many authors, notably 
by Mainardi et al. (2001, 2005) for the fundamental solution of the space-time 
fractional diffusion equation, and Saxena et al. (2006a, b) for fractional reaction- 
diffusion equations. The advantage of using Riesz-Feller derivative lies in the 
fact that the solution of the fractional reaction-diffusion equation containing this 
derivative includes the fundamental solution for space-time fractional diffusion, 
which itself is a generalization of neutral fractional diffusion, space-fractional 
diffusion, and time-fractional diffusion. These specialized types of diffusion can 
be interpreted as spatial probability density functions evolving in time and are 
expressible in terms of the H-functions in compact form. 



1 Introduction 

The review of the theory and applications of reaction-diffusion systems is con- 
tained in many books and articles. In recent work authors have demonstrated 
the depth of mathematics and related physical issues of reaction-diffusion equa- 
tions such as nonlinear phenomena, stationary and spatio-temporal dissipative 
pattern formation, oscillations, waves etc. (Frank, 2005; Grafiychuk, Datsko, 



f 



and Meleshko, 2006, 20076). In recent time, interest in fractional reaction- 
diffusion equations has increased because the equation exhibits self-organization 
phenomena and introduces a new parameter, the fractional index, into the equa- 
tion. Additionally, the analysis of fractional reaction-diffusion equations is of 
great interest from the analytical and numerical point of view. 

The objective of this paper is to derive the solution of an unified model of 
reaction-diffusion system (14), associated with the Caputo derivative and the 
Riesz-Feller derivative. This new model provides the extension of the models 
discussed earlier by Mainardi, Luchko, and Pagnini (2001), Mainardi, Pagnini, 
and Saxena (2005), and Saxena, Mathai, and Haubold (2006a). The present 
study is in continuation of our earlier work, Haubold and Mathai (1995, 2000) 
and Saxena, Mathai, and Haubold (2006a, 2006b). 



2 Results Required in the Sequel 



In view of the results 



J-i/2{x) = \ — cosx. 
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and (Mathai and Saxena, 1978, p. 49), the cosine transform of the H-function 
is given by 
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where R e [p + ^ m (h-)} > 0,Re[p + ^< n (^fj] < 0, \arga\ < ±7rQ,Q 
0; 

k > and Q = J2T=i B 3 E?= m+ i B 3 + £ "=i A j £^ = „ +1 Aj. 

The Riemann-Liouville fractional integral of order v is defined by (Miller and 

Ross, 1993, p. 45; Kilbas et al., 2006) 



,DrN(x,t) = f ^j\t-ur- 1 N(x, 



u)du, 
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where Re(v) > 0. 

The following fractional derivative of order a > is introduced by Caputo 
(1969; see also Kilbas et al., 2006) in the form 
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where §^f(x,t) is the m th partial derivative of f(x,t) with respect to t. 

The Laplace transform of the Caputo derivative is given by Caputo (1969; 
see also Kilbas et al., 2006) in the form 

m—l 

L { D?f(x, *);«} = s a F(x, s) - ^ s a - r -\f (r \x, 0+), (m - 1< a < m). (6) 

Following Feller (1952, 1971), it is conventional to define the Riesz-Feller 
space-fractional derivative of order a and skewness 9 in terms of its Fourier 
transform as 

F{ x D%f(x);k} = -* e a (k)f*(k), (7) 

where 

^ 9 a (k) = \k\ a exp[i(signk)—}, < a < 2, \6\ < min{a,2- a} . (8) 

When 6 = 0, then (8) reduces to 

F{ x Dgf(x);k} = -\k\ a , (9) 
which is the Fourier transform of the Weyl fractional operator, defined by 

D* m= 1 dU f* ^ du (10) 

This shows that the Riesz-Feller operator may be regarded as a generalization 
of the Weyl operator. 

Further, when 9 = 0, we have a symmetric operator with respect to x that 
can be interpreted as 

This can be formally deduced by writing — (k) a = — (fc 2 ) Q / 2 . For < a < 2 and 
|0| < min{a, 2 — a}, the Riesz-Feller derivative can be shown to possess the 
following integral representation in the x domain: 



+ sin[{«-6)K/2]j™ f{x ® +a f{x) dd}. (12) 

Finally, we need the following property of the H-function (Mathai and Sax- 
cna, 1978) 
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3 Unified Fractional React ion- Diffusion Equa- 
tion 



In this section, we will investigate the solution of the reaction-diffusion equation 

(14) under the initial conditions (15). The result is given in the form of the 
following 

Theorem. Consider the unified fractional reaction-diffusion model 

D?N(x, t) = Vx D%N(x, t) + *{x, t), (14) 

where rj, t > 0, x G r; a, 9, (3 are real parameters with the constraints 
< a < 2, \9\ < min(a, 2 — a), < (3 < 2, and the initial conditions 

N(x,0) = f(x),N t (x,0)=g(x));hv x G R, l $_+ ±a3 N(x, t) = 0, t > 0. (15) 

Here N t (x,0) means the first partial derivative of N(x,t) with respect to t 
evaluated at t = 0, 77 is a diffusion constant and $(x,i) is a nonlinear function 
belonging to the area of reaction-diffusion. Further x Dg is the Riesz- Feller 
space-fractional derivative of order a and asymmetry 8. o^f is the Caputo 
time-fractional derivative of order (3. Then for the solution of (14), subject to 
the above constraints, there holds the formula 

N(x,t) = — / f*(k)E /3A (-r ] t^ e a (k))exp(-ikx)dk (16) 

1 r°° 

+ 7T- / tg*{k)E p2 {-T 1 k a t^ e a (k))exp{-ikx)dk 
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+ 7T / / ^*(k,t-()E^ l3 {-r 1 k a t^ e a (k))exp(-tkx)dk. 

2lT JO J -00 

In equation (16) and the following, E a jj(z) denotes the generalized Mittag- 
LcfHer function (Saxena, Mathai, and Haubold, 2004; Berberan-Santos, 2005; 
Chamati and Tonchev, 2006). 

Proof. If we apply the Laplace transform with respect to the time variable t, 
Fourier transform with respect to space variable x, and use the initial conditions 

(15) and the formula (7), then the given equation transforms into the form 

s p N~(k,s) -s^V*(fc) -sP- 2 g*(k) = -r/*^(fc)iV~(fc,s) + $~(fc,s), 

where according to the conventions followed , the symbol ~ will stand for the 
Laplace transform with respect to time variable t and * represents the Fourier 
transform with respect to space variable x. 
Solving for N~ , it yields 

( ' S) ~ «0+r,tt«(jfe) + s0 + r,^ a (k) + sP + ^Hky [ ' 
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On taking the inverse Laplace transform of (17) and applying the formula 

L ^ l^T^} = t a -PE a>a _ p+1 {-at a ), (18) 

where Re(s) > 0, Re(a) > 0, Re(a — (i) > —1; it is seen that 

N*(k,t) = f%k)E^-Vt fj <(k))+g*(k)tE^ 2 (-Vt f3 <(k)) 

Jo 

The required solution (16) is now obtained by taking the inverse Fourier trans- 
form of (19). This completes the proof of the theorem. 



4 Special Cases 

When g(x) = 0, then by the application of the convolution theorem of the 
Fourier transform to the solution (16) of the theorem, it readily yields 
Corollary 1. The solution of the fractional reaction-diffusion equation 

g^N{x, t) ~ Vq^N(x, t) = t), x e R, t > 0, r, > 0, (20) 
with initial conditions 

N(x, 0) = f(x), N t (x, 0) = for x e R, 1< /3 < 2,£™ ±<J0 N(x, t) = 0, (21) 

where r] is a diffusion constant and &(x,t) is a nonlinear function belonging to 
the area of reaction-diffusion, is given by 

N(x,t) = f Gi{x-T,t)f{r)dT 
Jo 
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and 



1 

G 2 (x,t) = — / exp{-ikx)E fi4j {-f 1 t^ 9 a {k))dk 

, (a > 0). (24) 



alxr 3 ' 3 



(l,l/«),(/9,/3/«),(l,p) 
(l,l/a),(l.l).(l.P) 



In deriving the above results, we have used the inverse Fourier transform formula 

where Re(0) > 0, Re(-f) > 0, which can be established by following a procedure 
similar to that employed by Mainardi, Luchko, and Pagnini (2001). Next , if 
we set f(x) — S(x),$ — 0,g(x) — 0, where S(x) is the Dirac delta-function, 
then we arrive at the following interesting result given by Mainardi, Pagnini, 
and Saxena (2005). 

Corollary 2. Consider the following space-time fractional diffusion model 

=VxD?N(x, *),r?>0, x eR, 0</3<2, (26) 

with the initial conditions N(x,t = 0) = S(x),N t (x,0) = 0%™ ±ca N{x,t) = 
where 77 is a diffusion constant and S(x) is the Dirac delta-function. Then for 
the fundamental solution of (26) with initial conditions, there holds the formula 

N ( x t) - 1 gVr M |(l.l/°).(l,/Va),(M 1 (27 \ 

where p = 

Some interesting special cases of (26) are enumerated below. 

(i) We note that for a = (3, Mainardi, Pagnini, and Saxena (2005) have 
shown that the corresponding solution of (26), denoted by N%, which we call as 
the neutral fractional diffusion, can be expressed in terms of elementary function 
and can be defined for x > as 

Neutral fractional diffusion: 0<a~/3<2;9< min {a, 2 — a} , 

N o (x) = i^lWM (28) 

a{ ' 7rl + 2x a cos[(7r/2)(a-6)]+x 2a ' { ' 

The neutral fractional diffusion is not studied at length in the literature. 

Next we derive some stable densities in terms of the H-functions as special 
cases of the solution of the equation (26) 

(ii) If we set = 1, < a < 2; 9 < min {a, 2 — ajthen (26) reduces to space 
fractional diffusion equation, which we denote by L 9 a {x) is the fundamental 
solution of the following space-time fractional diffusion model: 

dN £' ~ =V x D%N(x, t), n>0,xeR, (29) 
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with the initial conditions N(x, t = 0) = 5(x), l *™ ±oa N(x, t) = 0,, where 77 is a 
diffusion constant and 5(x) is the Dirac-delta function. Hence for the solution 
of (29) there holds the formula 



L e a {x) 



1.1 



a 



(^t)l/a H 2,2 



(fft) 
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, < a < 1,|0| < a, 



(30) 



where p = The density represented by the above expression is known as 

a-stable Levy density. Another form of this density is given by 



1 



-H. 



2,2 



(Tit) 1 '' 



(l-i,4),(l-P,P) 
(0,1),(1-P,P) 



1< a < 2, \9\ < 2 -a, 



(31) 

(iii) Next, if we take a = 2,0 < (3 < 2,9 = 0, then we obtain the time 
fractional diffusion, which is governed by the following time fractional diffusion 
model: 

' A ,:,-./:.// •()../•, //.I)- i- 2, (32) 



dti 3 



with the initial conditions N(x,t = 0) = 6(x),N t (x,0) = 0^ ±oo N(x,t) = 
where 77 is a diffusion constant and 5(x) is the Dirac delta- function, whose 
fundamental solution is given by the equation 



(iv) Further, if we set a — 2, [3 — 1 and 
solution of the standard diffusion equation 



[ M 


(1,0/2) 


_(^)V2 


(1,1) 



(33) 

then for the fundamental 



§ i N( X ,t)=rj^N(x,t), 



with initial condition 



N(x,t = 0) = S(x), l ™ ±00 N{x,t)=Q, 



there holds the formula 



\ \x\ 


(1,1/2)" 


_^l/2 t l/2 


(1,1) 



(4 7 r^)- 1 /2 ea;p[ _^i_ ] 



(34) 



(35) 



(36) 



which is the classical Gaussian density. For further details of these special cases 
based on the Green function, one can refer to the paper by Mainardi, Luchko, 
and Pagnini (2001) and Mainardi, Pagnini, and Saxena (2005). 
Remark. Fractional order moments and the asymptotic expansion of the solu- 
tion (27) are discussed by Mainardi, Luchko, and Pagnini (2001). 
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Finally, for (3 = 1/2 in (14), we arrive at 
Corollary 3. Consider the following fractional reaction-diffusion model 

d] /2 N{x, t) = Vx D%N(x, t) + <f>(x, t), (37) 

where n, t > 0, x <G R; a, 9 are real parameters with the constraints 
< a < 2, \9\ < min(a, 2 — a), and the initial conditions 

N(x, 0) = f(x), for x e R, l ™ ±00 N(x, t) = 0. (38) 

Here 77 is a diffusion constant and <&{x,t) is a nonlinear function belonging to 

the area of reaction-diffusion. Further x Dg is the Riesz-Feller space fractional 

1/2 

derivative of order a and asymmetry 9 and D t is the Caputo time- fractional 
derivative of order 1/2. Then for the solution of (37), subject to the above 
constraints, there holds the formula 

N(x,t) = — J r{k)E 1/2 , 1 {-r 1 t^ e a {k))exp{-ikx)dk (39) 
1 r°° 

— C 1,2 di <5>*{kct-i)Ei_ i(-r ) k a t 1 - 2 ^ a (k))exp(-ikx)dk. 
2tt Jo J -00 2,2 



+ 



If we set 9 = in (39), then it reduces to the result recently obtained by the 
authors (2006a) for the fractional reaction-diffusion equation. 
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